Extremely-Metal Poor Stars in the Milky Way: 
A Second Generation Formed after Reionization 

Michele Trenti and J. Michael Shull 

University of Colorado, CAS A, Dept. of Astrophysical & Planetary Sciences, 389-UCB, 

Boulder, CO 80309, USA 

trentiOcolorado . edu, michael . shullOcolorado . edu 
ABSTRACT 

Cosmological simulations of Population III star formation suggest an initial 
mass function (IMF) biased toward very massive stars (M > 100 Mq) formed in 
minihalos at redshift ^ > 20, when the cooling is driven by molecular hydrogen. 
However, this result conflicts with observations of extremely-metal poor (EMP) 
stars in the Milky Way halo, whose r-process elemental abundances appear to be 
incompatible with those expected from very massive Population III progenitors. 
We propose a new solution to the problem in which the IMF of second-generation 
stars formed at 2; > 10, before reionization, is deficient in sub-solar mass stars, 
owing to the high cosmic microwave background temperature floor. The observed 
EMP stars are formed preferentially at 2; < 10 in pockets of gas enriched to 
metallicity Z > lO^'^'^ Zq by winds from Population II stars. Our cosmological 
simulations of dark matter halos like the Milky Way show that current samples of 
EMP stars can only constrain the IMF of late-time Population III stars, formed 
at 2; < 13 in halos with virial temperature T^i-^ ~ 10^ K. This suggests that pair 
instabihty supernovae were not produced primarily by this population. To begin 
probing the IMF of Population III stars formed at higher redshift will require 
a large survey, with at least 500 and probably several thousand EMP stars of 
metallicities Z ^ 10"^-^ Zp.. 



Subject headings: galaxies: high-redshift - Galaxy: evolution - intergalactic 
medium - stars: abundances - cosmology: theory - stars: formation 
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Introduction 



According to simulations, the first stars in the Universe began forming 50-100 miUion 
years after th e Big Bang (redshifts ^: = 30 — 50) in mini-halos with v i rial temperatures 



T ■ 



W K flTegmark et al.lll997l : lYoshida et al.ll2003l : iNaoz et al.l 120061 : iTrenti fc Stiavelli 



20071 ). To fix the semantics of early stellar populations, we define Population II (stars and 
gas) as material enriched to metallicities Z > Zcrit ~ at which point the star- 

formation mode and initial mass function (IMF) are altered through radiative cooling by 
metals, which dominate that of molecular hydrogen. The equality of the dynamical time 
and cooling time in a collapsing gas cloud marks the transition from a typical fragmentation 
mass of ~ 10^ Mq to a lower mass M < 10^ Mq, whose exact value d epends on possible 



2009 



couplin g with the Cosmic Microwave Background (CMB) temperature ( ISmith et al 
see also iBromm fc Loebl l2003l . lOmukai et al.ll2005l ). Therefore we define Population III as 



pristine, high-redsh i ft stars formed from gas at zero or near zero (Z < Zf^yit) meta llicity 



flBromm et al.ll200ll : iBromm fc LarsonI 12004 : ISantoro fc Shullll2006l : ISmith et al.l 120091 ) 



Initially very rare (one star per Gpc'^ at 2; ~ 50), these early generations of stars form 
at rates that increase rapidly in number toward l ower redshift until they are eventually sup- 
pressed by radiative and chemical feedback (see iTrenti et al.l |2009| and references therein) . 
The numerical simulations further suggest that the first-generation (Population III) stars 
formed in isolation, with an IMF with characteristic scale ~ 100 M^^ consistent with the 



fragmentation mass derived from simple cooling r nodels (see lAbel et al 



20031 : iBromm fc Larsonll2004j : lYoshida et al.ll2006l : IO'Shea fc Normanll2007L 120081 ). This char- 
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iBromm et al. 



acteristic mass is set primarily by the physics of hydrogen cooling, as these sta rs lack metals 
and dust. Molecular hydrogen (H2) can only cool the gas down to ~ 200 K (iGalli fc Palla 



19981 ). at which point the gas reache s a density ~ 10"^ cm ^ with corresponding Jeans mass 



see 



Mj ~ 10^ Mq (lYoshida et al.ll2006|). High acc retion rates onto the protostellar core (10 



O'Shea fc Norman 



20071 ) suggest that their fina l mass is likely to b e 



to 10-2 M0 yT-\ 

> 100 Mq. Radiative feedback will eventually shut down the accretion (iTan fc McKeell2004j ). 



but this phase has not yet been simulated self-consistently. 

In the presence of a strong flux in the H2 Lyman- Werner (LW) bands, with ener- 
gies between 11.15-13.6 eV, the primary (H2) coolant of min ihalos is photo-dissociated 
(ILepp fc Shulll Il983t iHaiman et al.l Il997l : iMachacek et al.l 120031 ) . Under these conditions, 
the minimum halo mass required for cooling increases from ~ 10^ Mc^ up to ~ 10^ M^t^ , 
reflecting a rise in vir ial temperature from Tvir ~ 10^ K to ~ 10^ K (IMachacek et al.ll2003l : 



Greif fc Brommll2006l ). In these larger halos, small clusters of Population III stars can form. 
The mass of the gas at the halo center is typically larger than in the case of a minihalo 
collapse, so a few tens of Jeans masses might be present. Because these second-generation 
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stars form relatively late, ai z < 15, some of them can be made out of reionized gas with 
correspondingly smaller masses (~ 40 Mq, see lYoshida et al.ll2006l ). 



A direct detection of Population 111 stars will be extremely challenging, even with the 
next generation of telescopes and instruments. A 100 Mq Population III star has an observed 
magnitude rriAB > 38 at ^ > 6, whereas the James Webb Space Telescope (JWST) is limited 
to rriAB ~ 31. Thus, JWST will only be able to observe such stars in the unlikely event 
they form in large clusters ( iTrenti et al.l 120091 ) . Supernovae from massive Population 111 
stars are bright enough to be seen at essentially any redshift, but their observed rate is low. 
Therefore, a combination of deep multi-epo ch imaging and large area coverage is needed 
dWeinmann fc Lillvll2005l : iTrenti et allboooh . 



Given the limited options for direct investigations of the properties of metal-free stars, 
indirect studies are the next best strategy. Such observations have been carried out to con- 
strain the IMF of Pop III stars by studying the properties of the most metal-poor stars in 



the lo cal universe ( iFreeman fc Bland-Hawthornll2002l : iBeers fc Christlieb 



20071 ). where remnants of the first Population III stars could be found (ITrenti et al. 



2005 



Frebe 



et al. 



20081). 



The idea behind these "galactic archeology" campaigns is that stars with very low metallic- 
ities Z < 10"^'^ Zq (approximately the critical metallicity marking the transition from 
Pop III to Pop II) are most likely second-generation objects formed from gas enriched 
by only one prev ious generat i on of stars. If this first generation of stars is very massive 



(M > 100 Mq), iHeger et al.l (|2003[ ) find that a significant fraction with initial masses of 



140 — 260 Mq will explode as pa i r instability supernova e (PISN) and produce no r-process 
elements (IHeger fc Woosleyll2002l : iTumlinson et al.ll2004l ). Surveys for extremely metal-poor 
stars in the Milky Way (MW) have found hu ndreds of stars with metallicity Z < 10~^Zq 



(IBeers fc Christhebl 12005 



Schorck et al.ll2009|). but so far none of those stars m atches the 



nucleosynthetic pattern of PISN progenitors (jTumlinsonll2006l : iFrebel et al.ll2009l ). Based on 
this observational evidence, it ha s been proposed that the IMF of Population III stars is 
peaked at lower masses, ~ 30 Mq (lTumlinsonll2006r ). This is in conflict with the results from 
numerical simulations of Pop III star-formation, which suggest higher masses for metal-free 
stars. 

Several past investigations studied the connection between observed m etallicity distri- 



bution of Galactic metal-poor stars and Population III stars. For example, iKarlsson et al. 
(120081 ) follow in detail the self-enrichment of a single ~ 10^ Mq halo and argue that second 
generation stars enriched by PISNe would reach a metallicity [Ca/H] ~ —2.5, higher than the 
cutoff used for current EMP searches. Hence, the lack of PISNe signatures in s tars considered 



to be second generation could simply be a result of an observational bias. iKarlsson et al. 



(120081 ) show that their results are insensitive to a redshift-independent change of the PopIII 
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star formation efficiency. However their work does not take into account the suppression of 
Population III star formation in minihalos due to LW feedback, which varies with redshift. 
Therefore the metalhcity overshoo ting by PISNe could p ossibly be less severe. A different 
explanation has been suggested by lSalvadori et al.l ( 120071 ). They use an analytic merger-tree 
code aimed at reproducing the z = Galactic metallicity distribution. As a result of fine tun- 
ing their free parameters to match the z = observations, they find a negligible probability of 
observing second generation stars. The robustness of their conclusion is difficult to evaluate, 
especially because their model has an idealized treatment of metal outflows and assumes a 
one-zone instantaneous mixing of metals for the gas outside star-forming halos. In addition. 
Population III stars are only assumed to form in halos with Tvir > 10^ K, complet e ly ne- 
glecting formation in mini halos cooled by if 2- Simila r conclusions to ISalvadori et al.l (120071 ) 
have also been reached by iKomiya et al.l (120071 . |2009| ) who also use an analytical merger-tree 
code (with instantaneous mixing approximation) and without taking into account the effects 
of radiative feedback on Population III star formation. 

Both of these physical effects are examined in this paper, with the goal of quantify- 
ing the IGM enrichment by Population III stars and predicting the expected fraction of 
EMP stars with PISNe signatures. The context of our study takes into account the com- 
plex interplay of radiative and chemical feedback during the reionization epoch. We employ 
high-resolution cosmological simulations of structure formation, both for an average region 
of the Universe (in a comoving volume of 10'^ Mpc^) and for three dark-matter halos sim- 
ilar in mass to the Milky Way. We post-process these simulations by taking into account 
additional physical pr ocesses tha t regu l ate star formatio i i, bey o nd those considered i n pas t 
investigations such as iTumlinsonI (120061 ) ; ISalvadori et al.l (120071 ) ; iKomiya et al.l (120071 . l2009l ) ; 
Karlsson et al.l (120081 ). These processes include: (1) radiative feedback on the cooling of 
metal-free minihalos; (2) spatial distribution of metal outflows from protogalaxies, missing 
in investigations based on analytical merger-tree codes; and (3) impact of the cosmic mi- 
crowave backgro und (CMB) temperature on the typical mass of Population II stars formed 



at high-redshift (IClarke fc Brommll2003l : iTumlinsonI 120071 : ISmith et al.ll2009h . 



This paper is organized as follows. In Section [2] we introduce our model and discuss 
the transition from Population III to Population II star formation and investigate its conse- 
quences for the chemical enrichment of second-generation gas in Section [31 Section [3] is based 
on analytical modeling and cosmological simulations representative of an average region of 
the high-redshift Universe. In Section H] we present results from three constrained-realization 
simulations of halos like the Milky Way, and we discuss the differences with the averaged 
results obtained for a random region of the Universe. Section \5\ concludes with a discussion 
of Milky Way EMP surveys needed to probe the metallicity distributions and set limits on 
PISN progenitors at z < 10. 
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2. Transition from Population III to Population II 



In this analy s is, we adopt the star formation model described in lTrenti fc Stiavellil (|2009[ ) 
and iTrenti et al.l (120091 ) . We consider the halo dynamics in a dark-matter scenario, based 
either on P ress-Schechter modeling plus an analytical expression for self-enrichment prob- 
ability (see iTrenti fc Stiavelli I2OO7I) or on detailed halo-merger histories derived from cos- 
mological simulations (ITrenti et al.l 120091 ) . Dark matter halos are populated according to 
an analytical cooling model that takes into account radiative feedback in the LW bands of 
H2, self-enrichment, and wind enrichment, the latter limited to the star formation histories 
from numerical simulations. A single, massive. Pop III star is assumed to form in halos 
with Z < Zcrit and lO^K ^ Tvir ^ lO^K. The precise value of Tvir required for Pop III for- 
mation depends on redshift and is derived considering the negative (H2 LW-band) feedback 
self-consistently determined from the sta r formation rate. Pop II I stars are assumed to end 
their life as a Pair Instability Supernova (iHeger fc Woosleyll2002[ ). while a constant fraction 
of gas (typically ~ 10"^) is converted into stars for metal-enriched halos with a minimum 
Tvir > 10^ K. Multiple Pop III stars are optionally allowed to form in the same halo when 
Lya cooling is possible, in metal-free halos with Tvir > 10^ K. 

Because of the large amount of meta l s (50-100 Mp-,) relea sed by a PISN explosion 
(iHeger fc WooslevI l2002f Ifleger et all boosl : ICal-Yam et"aDl2009h . a single Population III 
progenitor is sufficient to enr ich a descendant halo with Tvir ~ 10^ K to ^crit- The violent 
relaxation (ILynden-Belll Il967| ) associated with the virialization of the descendant halo pro- 
vides efficient mi xing of the metals, initially concentrated in a ~ 1 kpc region around the 
site of explosion (IBromm et al.ll2003[ ). However , the metal outflow from a PISN is conflned 
to the region occupied by a Tvir ^ 10^ K halo (IBromm et al.ll2003[ ). Hence, Population HI 
stars formed in isolation only contribute to genetic enrichment of their descendant halos 
and are not expected to pollute nearby halos. Metal outflows are instead expected from the 
Population II galaxies in halos with Tvir ^ 10^ K. These outflows are estimated to travel 
at average velocities of 40 — 100 km s~^ and can enrich a region of radius ~ 100 h^^ kpc 
( iMadau et al.l l2001t iFurlanetto fc Loebl 120031 : iTumlinson et al.l l2004j ) . We assume conserva- 
tively a star-formation efficiency = 0.01 and a metal yield of ym = 1/40, derived for a local 
z = IMF, hence likely an underestimate at 2; > 6 given the effect of the CMB temperature. 
We find that an outflow originating from a M = 10® Mq halo can enrich to Z ~ Zcrit a 
fully virialized 10® Mq halo within d < 33 kpc or a region containing the same mass 
in the process of virialization (with overdensity Ap/ (p) = 2) out to ci < 150 kpc. To 
model wind enrichment from the cosmological simulations, we therefore save snapshots with 
high-frequency {Az = 0.125 at 2; ^ 10 and Az = 0.25 at 10 ^ ^ ^ 15) and track during 
post-processing the expansion of metal outflows originating from Population II halos. We 
assume that a halo with Tvir ^ 10^ K is enriched to Z > Zcrit (if not already above this 
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threshold) once it is reached by the front of a metal outflow. 

Therefore, in our model, the gas is either pristi ne or above the critic al metallicity. For 



a lower value of Zcnt, possible in presence of dust ( ISchneider et al.ll2006l ). this assumption 
would be even better justified. Our standard analysis is carried out for Vwmd = 60 km s~^, but 
we also investigate Vwmd = 40 km s~^ and Vwind = 100 km s~^. Note that in our treatment 
of metal enrichment and star formation, metal-free halos with T^^,. ^ 10^ K are assumed 
to form a first generation of massive metal-free stars, immediately followed by a second 
generation of stars with Z > Zcrit- This implies that we are effectively setting the lifetime 
of a Population III star to zero and assuming instantaneous recycling for genetic enrichment 
in Tyir ^ 10^ K halos. In this approximation, we neglect the small fraction of halos that are 
quickly reached by a metal outflow, within a few Myr after hosting a burst of Population III 
star formation. Because winds propagate for > 100 Myr before reaching a metal-free halo, 
accounting for the lifetime of a massive Population III star would be equivalent to increasing 
Vwind by < 5%. This has a minimal impact in our model given the larger uncertainty on 
other parameters. 



We now summarize the main results of our model (ITrenti et al.ll2009l ). crucial to interpret 



the observed enrichment pattern in Galactic extremely metal-poor stars: 

1. Population III star formation in minihalos cooled by H2 starts at very high redshift 
{z > 45) for a comoving volume comparable to that of a MW-like halo. Radiative 
feedback begins quenching this mode of star-formation at z ~ 35 and completely 
suppresses it by z ~ 15. For a given star formation rate, the flux in the LW bands 
is almost IMF and metallicity independent, a very marked difference compared to the 
ionizing flux, which increases sign ificantly for a top-heavy IMF and for extremely low 



or zero metallicity (jSchaererl |2003| . Table 4). 



The metal-enriched star formation rate becomes higher than the metal-free rate at 
z ~ 25 and continues to increase almost exponentially until 2; ~ 10. At 2; ~ 20, Pop II 
stars are already the dominant star-formation mode and implies that metal-free stars 
are probably minor players in reionizing the Universe. At z < 25, LW photons are 
primarily produced by metal enriched stars. 

Population III star-formation continues to be possible at 2; < 15 at a comoving rate 
~ 10~^ M0 yr"'^ Mpc~^ in larger halos with Tyij. > 10^ K, insensitive to the LW photo- 
dissociating background. Winds from nearby Pop II halos enrich all the IGM by 2; ~ 4, 
progressively reducing the Pop III star formation rate, beginning at 2; ~ 10. 
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3. Second-generation Star Formation in average regions 



The complex star-formation history derived in iTrenti et al.l ( 120091 ) raises two questions 
in the context of the observed chemical abundances of EMP stars: First, what are the 
expected chemical progenitors of EMP stars? Are they Pop III stars in minihalos, Pop III 
stars in the larger Tyir ~ 10^ K halos (via self-enrichment), or Pop II stars (via winds)? 
Second, over which redshift range did the currently observable EMP stars form? 

We first address these questions for an average re gion of the Universe , using the = 
1024^ dark matter particle cosmological simulation of iTrenti et al.l ( l2009l ). with a comov- 
ing volume V = 10^ Mpc^, a single-particle mass 3.4 x 10"^ Mq, and a for ce resolution of 



0.16 /i ^ kpc. We adopt a cosmology based on the fifth-year WMAP data (IKomatsu et al. 
2009h : = 0.72, fi^ = 0.28, ilh = 0.0462, Hq = 70 km Mpc'\ and as = 0.817. 
The initial conditions have b een generated at z = 199 with a code based on the Grafic 



Eisenstein fc Hul ( 
friend halo finder ( 



algorithm f|Bertschingeii 120011 ) using a ACDM transfer function computed via the fit by 



9991) with spectral index = 0.96. Halos are identified with a friend-of- 



Davis et al. 



19851 ) using a linking length equal to 0.2 of the mean particle 



separation. A halo merger tree is constructed from the particle IDs to track metal enrichment 
derived from progenitor halos. 

We define "second-generation gas" as gas at Z > Zcnt, residing in a T^ir ^ 10"^ K halo 
with no progenitor that experienced Population II star formation in a previous snapshot. 
Next, we must distinguish between "second-generation gas" enriched by Pop II or Pop III 
stars. The earliest gas was enriched by Pop III stars in dark matter halos with T^^^ > 10^ K. 
During the halo merging history, such halos had at least one progenitor halo that hosted 
a Pop III star, but no progenitors with Pop II stars. In contrast, second-generation gas 
enriched by Pop II stars is defined as gas in a dark matter halo with Tvir > 10^ K whose 
progenitor halos had no star formation and were chemically enriched by winds originating 
in a nearby halc0. 

Figured] shows that most (about 60%) of the second-generation gas is enriched by Pop II 
winds. In addition. Pop III stars in minihalos (formed at z > 14) only contribute to a small 
fraction (about 10%) of second-generation enrichment; this gas is highlighted as the shaded 
area in Fig. [TJ This surprising result can be understood from the star-formation rate history 
plotted in Fig. 1 of ITrenti fc Stiavellil ( 120091 ). The Pop III star formation rate per unit time 
remains relatively constant from z 40 to z ^ 10, so that most of the stars are formed 



^The majority of halos with Tvir > lO** K in the simulation box have metal-enriched halos as progenitors 
and thus host at least third-generation stars. 
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at the lower end of the interval. (Note that cosmological time scales as t oc (1 + 

As time passes, metal-enriched winds from Pop II galaxies fill progressively more volume, 

enriching pockets of low-metallicity gas that never underwent self-enrichment because the 



stron g radiative (LW) background suppressed H2 cooling (see also Fig. 4 in iTrenti et al. 
2OO9I ). Figured] was obtained for fwind = 60 km s~^. Faster winds increase the fraction of 



the gas enriched by Pop II stars, while slow er winds favor a higher fraction of Pop III IGM 
enrichment (see Fig. 5 in lTrenti et al.ll2009l ). 



To translate the fraction of second-generation gas enriched by Pop III stars into a 
prediction for EMP star surveys, we need to address the second question raised above, 
namely when did currently observable EMP stars form? These stars can be observed only if 
they are bright nearby main-sequence stars or lie along the giant branch. Given that second- 
generation stars are quite old (they formed at 2; > 4 according to Fig. [1]), we can observe 
them only if their initial mass was M < 0.9 Mq. For the present-day IMF, most of the stars 
formed lie below this cut-off. However, this was not necessarily true at earl ier times, because 
coup l ing with the CMB radiation preven ts the gas to cool below Tcmb ( IClarke fc Bromm 
2o3; lTumlinsonl l2007l: ISmith et al.l l2009Hi_ 

In fact, the fragmentation scale of protostellar 
clouds is probably set by the Jeans mass ( lLarsonll2005l : iTumlinsonI 120071 ) . which depends on 
the temperature Tmin and density p of the gas: 

3/2 



Ml 



2mHG 



P 



-1/2 



'11 



By adopting a power-law equation of state, ITumlinsonI ( 120071 ) expresses the typical stellar 
mass in the IMF as a function of redshift: 

■max[2.73(l + z),8 



10 K 



(2) 



where a depends on the equation of state, with an expected range 1.7 < a < 3.35, and the 
normalization Mnorm is defined as: 

0.5 M© 



0.8" 



(3) 



based on a typical temp erature of 8 K and a typical mass of 0.5 M© observed in the local 
Universe ( |Kroupall2002l ). Equation |2 can t hen be used to define the IMF, parameterized in 
log-normal form (Eq. 4 in iTumlinson l2007[ l: 

2 



In 



dN 
d\nM 



A 



2a2 



In 



M 



M, 



c 



(4) 



^Note that the couphng is relevant only for gas at Z > Zait, because Population III halos cannot easily 
reach temperatures below T ~ 200 K, already higher than Tcmb at z < 70. 
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where A is a normalization constant and o" is a free parameter of order unity that controls 
the width of the distribution. 

From Eq. [2] and Eq. H] we can compute the number of second-generation stars (defined 
as stars formed out of "second-generation" gas) still observable today, with M < 0.9 Mq. If 
we are interested only in the relative fraction of EMP enriched by Pop III vs. those enriched 
by Pop II, our result does not depend on the EMP-star formation efficiency, as long as we 
assume a constant efficiency per unit mass in transforming gas to stars. Figure |2] shows this 
fraction, defined as 

e = ^Popiii —enriched —enriched + ^PopIII —enriched ) (5) 

as a function of {a, a). Here, a is the IMF mass width , and a is the exponent that regulates 
the typical stellar mass as a function of the CMB temperature (Eq. [2]). Even for a very 
broad IMF {a = 1.4) and a shallow redshift evolution of the typical stellar mass {a = 1.7), 
only up to ~ 25% of observable EMP stars are formed out of Pop Ill-enriched gas. For 
intermediate values of the IMF parameters, we obtain e ~ 0.15, while an extremely narrow 
IMF reduces e to ~ 0.06. These numbers imply that a significant sample of EMP stars with 
metallicity Z ~ is required to rule out the null hypothesis that Pop III stars are 

very massive and preferentially explode as PISN (see Tabled]). For a Poisson distribution, 
P{k) = N^^^exp {—Nexp)/k\, with expected value N^xp, the probability of zero events is 
P(0) = exp {—Nexp)- Thus, the minimum second-generation sample size to rule out that a 
fraction t] of Population III stars exploded as PISNe at a confidence level ^ (0 < ^ < 1) is: 

-log (1-0 ... 

For example, to rule out at 99% of confidence level (^ = 0.99) that 50% of Pop III stars 
explode as PISN, we need 41-135 second-generation stars, depending on the parameters 
assumed for the IMF of the EMP stars. 

With the same approach, we can estimate the size of the sample required to investigate 
the IMF of Pop III stars in minihalos at 2; > 15. The fraction of early-time Pop Ill-enriched 
second-generation stars for this case is shown in Fig. [31 We note the limited amount of 
gas enriched by this sub-class of Pop III stars (see Fig. [1]) as well as the increase of the 
typical EMP stellar mass as a function of redshift (see Eq. |2]) . It is highly unlikely to observe 
one of these second-generation stars today. In fact, we find e < 0.02 over the entire EMP- 
IMF parameter space. The main controlling parameter is a, the IMF width. Under the 
most optimistic assumptions for a and a in Table [H about 1775 second-generation stars are 
required to reject the hypothesis of 50% PISN from metal-free stars in minihalos at 99% 
CL. If the EMP-IMF is narrower, a sample size above 10^ sources is required. As we discuss 
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in 



Trenti et al.l (|2009[ ). the formation rate of minihalos, as measured in our cosmological 



simulation, might be shghtly underestimated owing to resolution issues. If this is the case, 
the size of the second-generation star sample is reduced accordingly. Conservatively assuming 
a factor-of-two higher Pop III formation rate in minihalos implies a sample size > 900 second- 
generation stars. 



4. Second-Generation Star Formation in Milky- Way like halos 

The star formation model in Section [3] describes an average region of the Universe at 
z > 4. Instead, the MW progenitors lived in a region with an above-average halo formation 
rate because of the presence of the large-scale overdensity that created the MW and the 
Local Group at z < 1. As a result, the star-formation history of the region is enhanced 
and tends to be shifted earlier in time compared with a typical region of the Universe. This 
introduces two competing factors that can affect the expected fraction of observable second 
generation stars enriched by Pop III stars. First, metal winds have had less time to reach 
nearby halos, so we expect a larger fraction of second generation gas enriched by Pop III 
stars. This effect is balanced by the higher CMB temperature floor present at earlier times, 
which makes it more difficult for stars formed out of this gas to have sub-solar masses and 
thus be observable today. 

To quantify this scenario and to obtain a detailed prediction to compare with EMP 
surveys in the MW, we consider a new set of cosmological simulations, with the same setup 
of the on e considered in Sec tion |3l except for the initial conditions. Here, we use the Grafic 



package (IBertschingerl 120011 ) and set the initial conditions to host a MW-like halo at z = 
at the center of the simulation box. Because of the prohibitive computational time needed 
to evolve to z = a. V = 10^ Mpc^ volume with = 1024^ particles, we first run a 
low-resolution [Nir = 128^) version of the simulation box to 2; = to measure directly the 
central halo mass realized in the initial conditions. We consider three different realizations 
of central halos to explore a range of masses and formation histories: Mi = 3.15 x 10^^ Mq 
(corresponding to 172, 981 low-resolution particles), M2 = 2.60 x 10^^ Mq, and M3 = 1.59 x 
10^2 Mq. From z = 1 to z = 0, halo 1 acquires 45% of its mass, similar to halo 3 (42%), 
while halo 2 experiences a major merger and acquires 59% of its mass. A projection of the 
halos at 2; = (region size ~ 1 Mpc^) from the low-resolution run is shown in the upper left 
panel of Figs. IHEl 

We then create the final full resolution (A^ = 1024^) box by keepi n g the same large-scale 



structure of the low-resolution version using the iHoffman &: RibakI ( 119911 ) method. Each 



low-resolution particle is thus associated with = 1024^/128^ = 512 "child" particles in 
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the high-resolution run. Based on the initial positions, we map the IDs of the particles in 
the central halo of the low resolution realization into the set of high-resolution particle IDs 
expected to be part of this central halo. We tested this method extensively with a series of 
constrained simulations at different resolution and found that, for a halo resolved with ~ 10^ 
low-resolution particles, the mass (and hence halo membership) is accurate at the > 97% 
level (Trenti et al., submitted). The setup we choose is s i milar in terms of mass resolution 
and numb er of partic l es to that adopted by Tumlinson ( 2009 lR with the only difference 



ution 



being that iTumlinsonI (120091 ) binned particle positions and velociti es to reduce the res o 
and complete the run from z = 4 to z = rather than using the iHoffman fc RibakI (Il99ll ) 
method. 

The post-processing analysis is the same discussed in Section [31 with the difference that 
the results on second-generation gas are given only for halos that are progenitors of the 
central z = MW-like halo. The comoving volume considered for normalization is that 
occupied at z = 199 by the particles that end up at 2; = in the central halo. The number 
of particles in this volume range from > 88 x 10^ for halo 1 to > 44 x 10^ for halo 3. 



Figures HHSl summarize the results for the different halos. The formation rate per unit 
redshift per unit volume of second generation gas is shown in the upper right panel of 
the figures as solid lines and compared to the box-averaged formation rate (dotted lines). 
As expected, the regions where the progenitors of MW-like halos live have an enhanced 
halo formation rate and higher clustering than an average region of the Universe. As a 
consequence, metal enrichment by Population III stars is more efficient (by > 50%) at 
z > 10. This is consistent with a very simple estimate based on linear theory: the region 
hosting a MW-like halo at 2; = has an average overdensity Ap/ (p) ~ 0.14 at z = 15. This 
in turn implies that the comoving n umber density of ha los with M = 1.4 x 10^ Mf?^ (hosting 



Popul ation III stars at z = 15, see iTrenti et al.ll2009l ). derived from the ISheth fc Tormen 



Jl999h mass function, is enhanced from ~ 5.6 Mpc ^ to ~ 8.9 Mp 



Because of the higher clustering, winds are also more efficient at polluting nearby, non 
self-enriched halos. At z < 10 the formation rate of Population III enriched gas becomes 
comparable to the average in the Universe. The enhanced formation rate of Tyir ~ 10^ K 
halos is balanced by a higher fraction of such halos reached by metal outflows. Consequently, 
the number density of wind-enriched halos also increases compared to the box average. By 
2; ~ 6 the number density of wind-enriched halos has returned to values similar to those of 
an average region of the universe. This can be again understood in terms of linear theory. 
Tvii- ~ 10^ halos are common at 2 < 6: based on their rarity, measured in terms of the Press- 



■^The preprint was posted on arxiv while this paper was being revised 
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Schechter variable v = S'^/a1)j^j{M), these halos have u ~ 1.38 on average in the Universe. 
The z = 5 overdensity associated with a MW-hke halo at z = only decreases this rarity 
to z/ ~ 1.14, and the number density of such overdensities is comparable: ~ 8.7 Mpc~^ on 
average versus ~ 9.4 Mpc~^ for the MW-like region. The three different halos we consider 
have all a similar enrichment history at z ^ 4. 

The predictions for the fraction e of second-generation stars enriched by Population III 
metals is shown in the bottom panels of Figures IHE] considering all Population III stars 
(left panels — equivalent to Figure and only Population III in minihalos (right panel 
— equivalent to Figure [3]). Given the two competing effects discussed above (higher halo 
formation rate and higher efficiency of wind pollution), the contour lines of e{a, a) are similar 
to those of the box-average. There is a moderate tendency to lower values for e, implying that 
the fraction of second-generation stars enriched only by Population III metals can be observed 
today is smaller than the one estimated for the average of the Universe. Tabled] contains the 
minimum number of stars required to rule out a 50% PISNe fraction in Population III stars 
for MW halo 1 based on Equation [61 Minimum sample sizes for halo 2 and halo 3 are similar. 
Depending on the IMF parameters, (a, a), the sample sizes for all Pop III enrichment range 
from a factor of a few difference with the box-average (at low a) to almost equal to the box 
average (for high a and low a). For enrichment by minihalos- Pop III in MW-like progenitors, 
the required sample sizes are about two times smaller than those derived in a random region 
(because of the factor two increase in the gas enriched by these stars compared to an average 
region of the Universe). 

Our main conclusion, namely the very small fraction of second-g eneration stars enriche d 



by Population III stars observable today, is similar to that found by lSalvadori et al.l ( 120071 ). 
although their modeling assumptions are different from ours. In particular, the instanta- 
neous, single-zone recycling of metals implies a sharp cut-off of Population III formation, 
while in our model we fin d a long tail of me tal-free star formation in halos with lower than 



average bias (see Fig. 3 in lTrenti et al.ll2009l ) 



Figure [7] explores the impact of Vu,ind on the results from our model for halo 1. The 
outflow velocity is the key parameter controlling the relative fraction of Population III to 
Population II enriched gas. In the absence of winds, only self-enrichment is possible, hence 
all second-generation gas would be enriched by Population III metals. On the other hand, 
faster propagation of metal outflows suppresses low-redshift Population III formation. For 
example, f^ind = 100 km s~^ reduces by > 75% the Population III star formation at z < 8 
compared to the standard Vwmd = 60 km s"^. Slower outflows, Vwind = 40 km s~^, increase 
the Population III star formation at low redshift by up to a factor two. This is reflected in 
the formation rate of second generation gas, shown in the upper panels of Figure. [3 The 
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lower panels contain the e{a, a) contour lines for these models. Naturally, the fraction of 
observable second-generation stars increases for lower wind speed, but not dramatically. For 
example, e < 0.23 in the standard model, e < 0.41 for Vwind = 40 km s~^, and e < 0.10 for 
Vwind = 100 km s-^ 



Conclusions and Discussion 



In this paper, we propose a novel solution to the conflict between the extremely top- 
heavy Pop III IMF predicted from cosmological hydrodynamic simulations and the lack of 
PISN signatures in the metal abundances of EMP stars observed in the Milky Way halo. 
Compared to past studies with a similar goal, our model takes into account additional 
physical ingredients relevant for star formation and chemical enrichment during the reion- 
ization epoch. The most important addition is the inclusion of radiative feedback in the UV 
Lyman- Werner bands of H2, bands, which quenches Population III star formation in mini- 
halo s and increases the min imum halo mass required for cooling to T^i,- ~10^Kby2;~15 
(see iTrenti &: Stiavellil |2009[ ) . Without this increase, we would have obtained a significantly 



higher self-enrichment efficiency and reduced the importance of metal outflows for enriching 
second-generation gas. We also base our analysis on high-resolution N-body simulations 
capable of resolving star forming minihalos, both in a 10^ Mpc^ (comoving) cosmological 
volume, as well as in three regions hosting a. z = halo similar in mass to that of the 
MW. These simulations allow us to quant ify metal outflows from proto-galaxies, improving 



the instantaneous one-zone assumpt ion in ISalvadori et al. 



merger trees, as well as the work of iKarlsson et al.l ( l2008l ). focused on self-enrichment of a 



(l2007l ). based on Press-Schechter 



single IO^Mq halo. A third novel key ingredient is the inclusion of the cosmic microwave 
background (CMB) coupling with cooling gas. This introduce s a floor in the cooling o f metal 
enriched halos and affects thei r initial stellar naass f unc tion (IClarke &: Bromnj|2003l ). This 
effect was also not modeled in ISalvadori et al.l ( 120071 ) or lKarlsson et al.l ( 120081 ). 



With our assumptions, we have shown that observations and numerical simulations can 
be consistent. In fact, we find that the majority of observable EMP stars with Z ~ 10~^'^ Zq 
were probably formed out of gas enriched to this critical metallicity by Pop II winds a.t z < 10. 
Pop III stars formed in minihalos cooled by H2 enrich only a small fraction, e ~ 10~^ to 
10~^, of the gas that ends up in EMP stars with M < 0.9 Mq. The fraction of observable 
second-generation stars formed in regions enriched by late-time Pop III stars is even higher, 
ranging from 6% to 25%, depending on the details of the modeling, in particular the IMF 
width of the EMP stars. 



The minimum sample sizes reported in Table [T] can be compared with the metallicity 
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distribution of EMP stars based on Galactic survey s such as t 



fIChristlieb et al.ll2008[ ) and the SDSS/SEGUE survey fllvezic et al. 



le Ha mburg/ESO survey 



2008 ) that have identified 



several hundred metal-poor stars with [Fe/H] < —3.0 (ISchorck et al.ll2009[ ). These surveys 
show a sharp cutoff at [Fe/H] < —3.6, consistent with the theoretical expectation of Z^^t ~ 
\^^^-^Zq. To infer the Pop III IMF, we should ideally use only those EMP stars close 
to the critical metallicity for the transition to Pop II star formation. This would help to 
avoid contamination in the sample by third-or-higher generation stars formed in gas with 
a lower than average enrichment history. This issue is especially important if the lack of 
PISN chemical signatures in second-generation stars is used to argue against very massive 
Pop III star progenitors. A single PISN producing 10 Mq of iron can pollute 10^ Mq of 
gas (the typical content of a T^r ~ 10^ K halo) to a metallicity [Fe/H] ~ —3.2, assuming 
homogeneous mixing. This is a level slightly higher than than the observed cut-off at [Fe/H] 
~ —3.6, suggesting that it would be important to improve instantaneous one-zone metal 
enrichment assumptions. Except for this caveat, we thus suggest that the analysis should be 
confined to a sample of EMP stars w ith [Fe/H] < —3.4 t o identify second-generation stars. 
The HES survey has 14 such stars (ISchorck et al.l 12009). A few oth er EMP with [Fe/H] 
< —3.5 also have published metallicities (see iBeers fc Christliebl l2005l for a recent review). 
Based on this sample size, we would have expected on average a few of these stars to show 
PISN signatures in their abundances if a significant fraction of late-time Pop HI stars were in 
the 140-260 M0 range. A larger sample of EMP stars, likely available in the near future with 
extensions and/or follow-ups of the SEGUE and HES surveys, will be able to set stronger 
limits on the fraction of PISN originating at z < 10. Investigating the shape of the IMF for 
Pop HI stars in minihalos would require an increase of about two orders of magnitude in the 
current sample of EMP stars. 

Since the majority of Pop HI numerical investigations have focused on Pop HI stars in 
minihalos, there is no inconsistency with the observed abundances of EMP stars. The sample 
size used to search for PISN signatures is simply too small. The current EMP sample is now 
approaching a size that can yield solid constraints on the IMF of late-time Pop HI stars (see 
Tabled]). The absence of PISN signatures is still lacking a significant confidence level, but it 
is consistent with the predictions for Pop HI st ar formation in (pa. rtially) reionized gas. It 
points toward a typical IMF mass of ~ 30 M© ( Yoshida et al. 20061 ). 



Our scenario for the formation of second-generation stars is also consistent with th e 
observed distribution of the [Mg/H] ratio for EMP stars (see Fig. 1 in iFrebel et al.l l2009l ) . 
Magnesium abundances below the rninimum level predi cted by PISN nucleosynthesis, [Mg/H] 
>, —6.1 see iHeger fc Wooslevll2002l : iFrebel et al.ll2009l ) are naturally expected if the enrich- 



ment of second-generation stars with M < 140 Mq is driven by Pop II winds and/or by 
Pop HI stars formed in (partially) reionized gas. 
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From our Figures illustrating the chemical enrichment of MW-like halos, we see 
that the comoving star formation rate per unit redshift averages several times 10^ Mq Mpc~^ 
over a redshift interval Az ~ 6. Thus, in a comoving volume of 20 Mpc^, approximately 
equivalent to the region that collapses into a MW-like halo, we estimate that ~ 2.7 x 10^ M© 
of second-generation gas was available to form stars. This is ~ 2% of the total baryonic 
mass of the Milky Way today. We expect that several thousand second-generation EMP 
stars could be observable today. For example, adopting a rather inefficient star formation 
efficiency of 0.1% and assuming that only 10~^ of the second-generation stars had an initial 
mass below 0.9 Mq, we still expect a total of 2 x 10^ EMP stars with [Fe/H] ~ —3.5 in the 
MW halo. 

Finally, we have shown that the transition from Population III to Population II star 
formation that occurred at 2; > 4 is not too different in the region that collapses at z = 
in a MW-like halo compared to the average in the Universe (Figures HUH]). At ^ > 10, 
the Population III star formation rate, and its associated chemical enrichment of second 
generation gas are about a factor two higher compared to the box-average. At z = 10, metal 
outflows from proto-galaxies are also more efficient, but by less than about a factor two. 
However, by 2; < 6, the formation rate of 10^ Mq halos has returned to the average value in 
the Universe. Our findings demonstrate that constrained realizations of MW-like halos are 
not necessarily required to investigate the first steps in the build-up of our own Galaxy. 

We thank Britton Smith and Massimo Stiavelli for useful comments and discussions. 
We acknowledge support from the University of Colorado Astrophysical Theory Program 
through grants from NASA (NNX07AG77G) and NSF (AST07-07474). 
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Fig. 1. — Formation rate per unit redshift of second-generation gas (M© Mpc~^ dz~^) en- 
riched to metallicity Z ~ 10~^'^ Zq as measured from a 10^ Mpc^ cosmological simulation 
that takes into account radiative feedback for Populat ion III formation, self-enrichment of 
halos, and metal winds propagating at 60 km s~^ (see iTrenti et al.ll2009[ ). The rate of gas 



enriched by Pop III stars is shown as a solid blue line, while the rate for Pop Il-enriched gas 
is shown as solid red line. The blue shaded area represents the gas enriched by Population 
III stars formed in minihalos (cooled by II2). Most second-generation stars in the volume 
are formed near the end of the reionization era, at redshifts 2; < 10. 
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Fig. 2. — Upper panel: Image showing the ratio e of stars enriched by Population III to total 
number of second-generation stars (with initial mass M < 0.9 M©) on the main sequence or 
giant branch in the current universe, plotted as function of log-normal IMF parameters a 
and a. Color coding for e is given by the color bar on the right. Lower panel: contour lines 
for e{a, a) = 0.07, 0.1, 0.15, 0.20, 0.25, 0.29 (from bottom to top). 
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Fig. 3. — Same as Fig. [21 but considering only second-generation stars enriched by the 
Population III stars in minihalos. This second-generation gas is represented as shaded area 
in Fig. [H Contour lines are for much smaller values of e = 10~^, 10~^, 10~^, 10~^, 2 x 10~^ 
from bottom to top. Studies of extremely metal-poor stars will need a very large sample to 
constrain the IMF of Population III stars progenitors in minihalos. 
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Fig. 4.— MW-like halo 1 (Mi = 3.15 x 10^^ Mq at z = 0). Upper left: halo projection 
at z = (from the low-resolution run). The image size is 1.3 x 0.97 Mpc^. Upper right: 
formation rate of second-generation gas at high redshift within the comoving volume that 
ends up in the halo at 2; = (solid lines). This panel is the equivalent of Fig. [H obtained 
for Vy^ind = 60 km s~^. For comparison, dotted lines show the average formation rate of 
second-generation gas from Fig. [H Smoothing over two data points has been applied for the 
solid lines. Bottom left panel: Contour lines for the ratio e of stars enriched by Population 
III to total number of second-generation stars for halo 1. This panel is the equivalent of 
Fig. 121 Bottom right panel: like in bottom left panel, but considering only Population III 
stars in minihalos (equivalent to Fig. [3]). 




Fig. 5.— Same as Figure S for MW-like halo 2 (Ms = 2.60 x 10^^ M© at ^ = 0). 




Fig. 6.— Same as Figure S for MW-like halo 3 (M2 = 1.59 x 10^^ M© at ^ = 0). 
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Fig. 7. — Influence of metal outflow speed on our modeling for MW-like halo 1 (see Fig- 
ure H]). Left column shows results with Vwind = 100 km s^^, and right column shows 
Vwind = 40 km s~^. Upper panels (formation i rate of second-generation gas) show for com- 
parison our standard model, f^md = 60 km s~^, as dashed lines. Bottom panels show contour 
lines for the ratio e of stars enriched by Population III to total number of second-generation 
stars. 
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Table 1: Required sample size of EMP stars with Z ^ 10 ^-^ 



a 


a 


^min— 


AllPopIII 


^min— AlinihaloPopI 1 1 


1.7 


0.3 


118 


274 


> W 


> 10^ 




1.7 


0.6 


65 


95 


1.35 X 10^ 


6.63 X 


10^ 


1.7 


1.0 


41 


49 


1775 


1008 




2.5 


0.3 


130 


449 


> 10^ 


> 10^ 




2.5 


0.6 


99 


202 


> 10^ 


> 10^ 




2.5 


1.0 


62 


87 


5.09 X 10^ 


2.54 X 


10^ 


3.35 


0.3 


135 


588 


> 10^ 


> 10^ 




3.35 


0.6 


120 


344 


> 10^ 


> 10^ 




3.35 


1.0 


86 


151 


4.61 X 10^ 


2.30 X 


10^ 



Note. — Minimum nmxibcr {Nmin, see Equation [5]) of second-generation stars needed to rule out at a 
confidence level of 99% the hypothesis that a fraction 50% of Population III stars exploded as PISN (third 
column) as a function of the second-generation IMF parameters a and a. The first value of Nmin refers 
to the box average, the second value to MW-like halo 1. The last column reports the minimum number of 
second-generation stars needed to rule out, at the same confidence level, PISN in minihalos Population III 
stars. For some of the IMF parameter space, the number of s tars needed excee ds 10''. The values of the 
characteristic IMF width a reported here are those assumed in Tumlinson ( 20061 ). 



